Chiral symmetry patterns of excited mesons with the Coulomb-like linear 

confinement. 
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The spectrum of qq mesons in a model where the only interaction is a linear Coulomb-like instan- 
taneous confining potential is studied. The single-quark Green function as well as the dynamical 
chiral symmetry breaking are obtained from the Schwinger-Dyson (gap) equation. Given the dressed 
quark propagator, a complete spectrum of "usual" mesons is obtained from the Bethe-Salpeter equa- 
tion. The spectrum exhibits restoration of chiral and U(1)a symmetries at large spins and/or radial 
excitations. This property is demonstrated both analytically and numerically. At large spins and/or 
radial excitations higher degree of degeneracy is observed, namely all states with the given spin fall 
into reducible representations [(0,1/2) (1/2,0)] x [(0,1/2) (1/2,0)] that combine all possible 
chiral multiplets with the given J and n. The structure of the meson wave functions as well as the 
form of the angular and radial Regge trajectories are investigated. 

PACS numbers: 11.30.Rd, 12.38.Aw, 14.40.-n 



INTRODUCTION 

In spite of numerous efforts for 30 years the physics 
which drives the structure of hadrons in the light fla- 
vor sector is not yet completely understood. It is clear, 
however, that here the most important elements of QCD 
are spontaneous breaking of chiral symmetry and con- 
finement. The hadron spectrum is a key to understand 
QCD in the nonperturbative regime, the role and inter- 
relations of confinement and chiral symmetry breaking. 

The spontaneous breaking of chiral symmetry is cru- 
cial for a proper description of the physics of the low- 
lying hadrons. It is obvious, for example, in the case of 
a pion, which has a dual nature. From the chiral sym- 
metry point of view it is a (pseudo) Goldstone excitation 
associated with the spontaneously broken axial part of 
chiral symmetry [l[. From the microscopical point of 
view, it is a highly collective quark-antiquark mode Q. 
The latter feature persists in any known microscopical 
approach to chiral symmetry breaking that is consistent 
with the conservation of the axial vector current. By now 
it is well understood that both these views are comple- 
mentary and consistent with each other. That the chiral 
symmetry breaking is crucially important for a proper 
physics of the nucleon follows e.g. from the large pion- 
nucleon coupling constant or from the Ioffe formula Q 
that relates, though not rigorously, the nucleon mass to 
the quark condensate of the vacuum. 

It has been realized in recent years that the physics in 
the upper part of the light hadron spectra is quite dif- 
ferent. In this case the chiral symmetry breaking in the 
vacuum is almost irrelevant, as evidenced by the persis- 
tence of the approximate multiplets of chiral and U(1)a 
groups both in baryon 0J1, [i| and mesons 0] spectra, 
for a short overview see [8(. This phenomenon, if con- 
firmed experimentally by discovery of still missing states, 
is referred to as effective chiral symmetry restoration or 



chiral symmetry restoration of the second kind. It should 
not be mixed up with the chiral symmetry restoration at 
high temperatures and/or densities. In the latter case the 
system undergoes a phase transition and the chiral order 
parameter vanishes - the quark vacuum becomes trivial. 
In contrast, nothing happens of course with the quark 
condensate of the vacuum once we are interested in the 
given isolated hadron at zero temperature. If a hadron 
is highly excited, this quark condensate of the vacuum 
becomes simply unimportant for physics, because the va- 
lence quarks decouple from the condensate. 

A fundamental origin of this phenomenon is that in 
highly excited hadrons the semiclassical regime is mani- 
fest and semiclassically the quantum fluctuations of the 
quark fields are suppressed relative to the classical contri- 
butions that preserve both chiral and U(1)a symmetries 
Q . This general claim has been illustrated [§] within the 
exactly solvable chirally symmetric confining model [l(| , 
which can be considered as a generalization of the large 
N c 't Hooft model in 1+1 dimensions ll| to 3+1 dimen- 
sions. The key issue is that the chiral symmetry break- 
ing Lorentz-scalar dynamical mass of quarks arises via 
loop dressing and vanishes at large momenta, where only 
classical contributions survive. When one increases the 
excitation energy of a hadron, one also increases the typ- 
ical momentum of the valence quarks. Consequently, the 
chiral symmetry violating dynamical mass of quarks be- 
comes small and chiral symmetry gets approximately re- 
stored in the given hadron. This theoretical expectation 
has been proven by the direct calculation of the heavy- 
light excited meson spectra with the quadratic confining 
potential and of the light-light meson spectra with 
the linear potential 

The effective chiral restoration has transparently been 
illustrated within a simple toy model at the hadron level 
that generalizes the sigma model and contains an infinite 
amount of excited pions and sigmas 14( . 
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FIG. 1: Graphical representation of the Schwinger-Dyson 
(A) and Bethe-Salpeter equations (B) in the ladder approxi- 
mation. 



Restoration of chiral symmetry requires that hadrons 



should decouple from the Goldstone bosons 
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A hint for such a decoupling is well seen: 
The coupling constant for the decay process h* — > h + 
7r decreases very fast once one increases the excitation 
energy of a hadron. 

The effective restoration of chiral symmetry has been 
discussed recently within different approaches in refs. 

In this paper we present a formalism that leads to the 
meson spectrum reported in our recent Letter [l3| as 
well as some additional numerical results. We also study 
meson wave functions and clarify origins of a very fast 
restoration with increasing J and of a slow restoration 
with increasing the radial quantum number n. 

Our choice of a model is motivated by the following 
constraints. The model must be (i) 3+1 dimensional, 
because we know that in 1+1 dimensions the effective 
chiral restoration does not occur due to absence of the 
rotational motion and spin; (ii) chirally symmetric; (hi) 
field-theoretical in nature in order to be able to exhibit 
spontaneous breaking of chiral symmetry; (iv) confining; 
(v) preserving the conservation of the axial vector cur- 
rent; and (vi) solvable. Then, upon solving a spectrum 
within such a model one can judge whether the chiral 
symmetry patterns appear or not and whether the phys- 
ical arguments presented above are correct. 

The only known model that satisfies all these con- 
straints is the Generalized Nambu and Jona-Lasinio 
model of ref. [loj ]. In this model the only interaction 
between quarks is the instantaneous confining Lorentz- 
vector potential between color quark charges. This model 
belongs intrinsically to the class of large N c models, be- 
cause only ladders are taken into account upon solving 
the Schwinger-Dyson and Bethe-Salpeter equations, see 
Fig. 1, and there are no vacuum fermion loops. Differ- 
ent aspects of this model, in particular chiral symmetry 
breakin g, h ave been studied in many works, see for exam- 
ple mr 
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28|| and the spectrum of the lowest exci- 
tations of the light-light mesons with the linear potential 
has been calculated in ref. [29| . It has been recently ap- 
plied to study confinement in the diquark systems [301 ] . 



From the underlying physics point of view, this model 
can be directly related to the Gribov-Zwanziger scenario 
of confinement in Coulomb gauge 31]. In Coulomb gauge 
a Hamiltonian setup to QCD can be used and nu- 

merous studies in different approaches suggest that the 
Coulomb interaction in the infrared turns out to be a con- 
fining linear potential, see e.g. 

ESS, Hit Actually 
the infinitely increasing Coulomb potential is a neces- 
sary but not a sufficient condition for confinement mech- 
anism in QCD and the strength of the Coulomb potential 
must be larger or equal to that extracted from the Wilson 
loop 34( . Then there is no contradiction with the recent 
observation on the lattice that in the deconfining phase 
above T c a strong Coulomb confinement still persists [38| . 
This shows, however, that the actual confining scenario 
in QCD is richer than a simple linear Coulomb-like po- 
tential. Hence some elements of QCD are still missing 
in this picture and one needs them in order to obtain a 
completely realistic picture of hadrons. 

However, our purpose is different. We use this model to 
answer a principal question whether effective restoration 
of chiral symmetry in excited hadrons occurs or not and 
if it does - to use the model as a laboratory to get an 
insight. This model satisfies all the criteria (i)-(vi) and 
hence is completely adequate for this purpose. 



THE HAMILTONIAN AND THE INFRARED 
REGULARIZATION 

The GNJL model is described by the Hamiltonian [l(| 

+ \ J d 3 xd 3 y J°{x,t)K? v (x- y)4(y,t), (1) 
with the quark current-current interaction 



A" 



(2) 



parametrized by an instantaneous confining kernel 
K?~ (x — y) of a generic form. We use the linear con- 
fining potential, 

Kf v (x-v) = g^g^vdi-gi), (3) 

and absorb the color Casimir factor into the Coulomb 
string tension, 



A a A Q 



V(r) 



(4) 



The Fourier transform as well as the loop integrals cal- 
culated with the linear potential are infrared divergent. 
Hence we have to perform the infrared regularization and 
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suppress the small momentum behavior of the linear po- 
tential by introducing a cutoff parameter into the poten- 
tial. Then in the final answer this cutoff parameter must 
be sent to 0, i.e. the infrared limit must be taken. While 
quantities, such as the single quark Green function, can 
be divergent in the infrared limit, which means that a sin- 
gle quark cannot be observed, we have to check that any 
color-singlet observable, like meson mass, is finite and be 
not dependent on the infrared regulator in the infrared 
limit. 

There are slightly different and physically equivalent 
ways to perform the infrared regularization in the litera- 
ture. We define the potential in momentum space as in 
ref [H 



V(p) = 



87TCT 



(p 2 + Mm) 2 ' 



(5) 



Then, upon the transformation back into a configura- 
tional space 
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(2tt) s 



Spe l P p V(p) = - 



crexp(-/ii R r) 



G 
MlR 



MIR 

+ ar + 0(fi m ) (6) 



one recovers that the potential contains the required lin- 
ear potential plus a term that diverges in the infrared 
limit nm — > 0. Note that the divergent part is canceled 
exactly in all physical observables and there is no depen- 
dence on the particular way of the infrared regularization. 

Given this prescription all the integrals are well de- 
fined and we can proceed by solving the Schwinger-Dyson 
(gap) and Bethe-Salpeter equations. Note that there are 
no ultraviolet divergences if only a linear potential is 
used. They will persist, however, once a Coulomb po- 
tential is added. Then the ultraviolet regularization and 
renormalization would be required. Since our purpose is 
to study the chiral symmetry restoration in the highly 
excited light-light states, where a role of the Coulomb 
interaction is negligible, such a complication is not re- 
quired. 



THE GAP EQUATION AND THE CHIRAL 
SYMMETRY BREAKING 



Due to interaction with the gluon field the dressed 
Dirac operator for the quark becomes 

D(p ,p) = iS^ 1 (p ,p) = D {po,p) - £(p ,p), ( 7 ) 

where the bare Dirac operator with the bare quark mass 
to is 



and Tt(po,p) is the quark self energy. For an instanta- 
neous model the latter becomes energy (po) independent 
and is given by 



W) 4 



V{\p-qy)lo-^i 



S (q ,q)-'S{q) 



7o . (9) 



Inserting (J9j> into ([7|) yields an integral equation for 
S(po,p). With the ansatz 

Z(£) = A(p)-m + j-p(B( P )-p), (10) 

where here and in the following p means \p\ and p = p/p, 
the equation ([7]) can be cast into a coupled system of 
integral equations 



A { p) = m + \j^_V(k) 
B(p)= P+ \j^V { k)p 



M(q) 
Q(q) ' 



uj(q) 



(11) 
(12) 



where the Lorentz-scalar A(p) and Lorentz-spatial vector 
B(p) parts of self-energy contain classical and quantum 
loop contributions [§]. Here 



M(p)=p 



B(p) 



(13) 



is called the quark mass function, uj(p) = \J M 2 (p) + p 2 
and k = p — q. The appearance of a dynamical mass 
signals the spontaneous chiral symmetry breaking. In the 
chiral limit to = it has exclusively a quantum (loop) 
origin. 

Notice that in (|11|12| the integrals are divergent at k = 
unless an infrared (IR) regulator /iir > is introduced. 
By defining the chiral angle (p p as 



sm< 



cos (p p 



A(p) 


M(p) 


uj(p) 


Q(p) 


B(p) 


P 


uj(p) 





where 



u J 2 (p) = A 2 (p) + B 2 (p), 



(14) 
(15) 

(16) 



(Illll2p can be transformed into one integral equation 

p sin ifp — m cos ip p = 

1 f A , s , 

2 J (2tt)3 ^ ' (smipqcosipp -p ■ qcosip q smip p ) 

(17) 

for the chiral angle. The integral in (jTTJ) is convergent at 
k = since the infrared divergences arising by integrat- 
ing over the two terms of the integrand mutually cancel. 
Hence the chiral angle and the dynamical mass 



Dq(pq,p) = iS 1 (p 0l p) = polo -p-i-m, (8) 



M(p) = pt&ntpp, 



(18) 



4 




This yields 



A(p) 
B(p) 



2mir 
a 

2/ilR 

a 
2mir 



sinipp + A f (p), 
cos ip p + Bf(p), 
•W/(P)> 



(22) 
(23) 
(24) 



FIG. 2: Mass function and chiral angle in the chiral limit for 
different values of the infrared regulator /iir. All quantities 
are given in appropriate units of \fa. 



where At(p), Bf(p), and are infrared- finite func- 

tions. 

Numerically the integration in (|17[) with = 
is feasible but demanding. Alternatively one can also 
solve (jlll 12[) with sufficiently small but finite values for 
//ir = 0. Still, in particular for very small values of /zir, 
special care has to be taken for the numerical integration 
in the vicinity of q = p. The convergence of M(p) for 
MiR ~~ * m t nc chiral limit (m = 0) has been demon- 
strated in Refs. [26|, [3(|. The mass function and the 
chiral angle, which are in agreement with the previous 
studies, are shown in the upper and lower plot of Fig. 1, 
respectively. The result for /xtr = 10~ 4 y/a is already a 
very good approximation of the infrared limit which can 
be determined numerically by extrapolating the results 
to /Uir = 0. 

An important issue is a momentum dependence of the 
chiral angle and dynamical mass. Their nonzero values 
are uniquely related in the chiral limit with the sponta- 
neous symmetry breaking. The effect of symmetry break- 
ing is large at low momenta and the chiral angle goes to 
7r/2 at zero momenta. At large momenta the mass func- 
tion and the chiral angle go to zero. This property is 
crucial for a proper understanding of chiral symmetry 
restoration in excited hadrons. 



are infrared finite. Then the quark condensate is given 
as 



(w) = T dp p 2 sin tp p . 

n Jo 



(19) 



On the other hand, A(p), B(p) and consequently also 
lu(p) diverge like fi^. The latter can be determined from 
A(j>) and B(p) with Eq. (JTHJ) or directly by combining 

cnn to 

uj(j>) = m sin ip p + p cos ip p + 

1 / d 3 q . , . . „ A . 

V (k) (sin ip q sm ip p + p ■ q cos ip q cos <p p ) . 



2 J (2tt) 3 

(20) 

The infrared divergences can be explicitly projected out 
by using 



wr^o 7T 2 j (k z + fj,f K y 

[ d\ 6(p-q)f(q) = f(p). 



(21) 



CHIRAL AND U(1) A SYMMETRY PROPERTIES 
OF THE NONINTERACTING TWO-QUARK 
AMPLITUDE 

Given a dressed quark Green function, one can obtain 
a spectrum of the quark-antiquark bound states from the 
homogeneous Bethe-Salpeter equation in the rest frame. 

Before solving the Bethe-Salpeter equation we would 
like to discuss a limiting case when the self-interaction 
is switched off, M(p) = 0. In such a case there is no 
self-energy and the quark condensate is identically zero. 
The chiral and U(1)a symmetries are not broken. Then 
the Bethe-Salpeter amplitude is reduced to a system of 
massless quarks with definite chirality. In such a case 
properties of the system are unambiguously determined 
by the Poincare invariance as well as chiral and U(1)a 
symmetries (tJ. The Poincare invariance prescribes the 
standard quantum numbers J PC . The chiral and U(1)a 
properties of the system then are uniquely determined by 
the index of a representation of SU(2)l x SU(2)r that 
is compatible with the given J PC and isospin /. Then a 
complete list of states is 
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(l/2,l/2) : 1,0-+ 
(l/2,l/2) 6 : 1,0++ 



0,0+ H 
0,0~ H 



(25) 



2k, k=l,2,. 



(0,0) 
(l/2,l/2)« 
(l/2,l/2) 6 
(0,1)0(1,0) 



0, J"" 

1, J- + 
1,J ++ 
1,J++ 



0, J++ 
0,J++ 

0, J~+ 

1, J 



(26) 



J = 2k-l, k=l,2,. 



(0,0) 


0, J++ <- 


-» 0, J 


(l/2,l/2) 


1,J + -- 


— * o,j— 


(l/2,l/2) 6 


1,J— ^ 


-o,j+- 


(0,1) ©(1,0) 


1,J— ^ 





(27) 



The sign < — ► indicates that both given states belong to 
the same representation and must be degenerate. Note, 
that the structure of the chiral multiplets is different for 
mesons with different spin. In particular, for some of the 
J PC the chiral symmetry requires a doubling of states. 

The U(1)a multiplets combine the opposite spatial 
parity states from the distinct (1/2, 1/2) Q and (1/2, l/2)& 
multiplets of SU(2) L x SU(2) R . 

The energy of the rotating massless quark is the same 
irrespective whether it is right or left, because there is 
no spin-orbit force in this case Q. Then one expects 
a degeneracy of all chiral multiplets with the same J. 
This means that the states with the same J fall into a 
reducible representation 



[(0,1/2)0(1/2,0)] x [(0,1/2)0(1/2,0)] 



(28) 



that combines all possible chiral representations of the 
quark-antiquark system with the same spin. 

Certainly, once the self-energy is switched on, the chi- 
ral symmetry gets broken and there cannot be exact chi- 
ral and U(1)a multiplets. The effective restoration of chi- 
ral symmetry in excited hadrons is defined to occur if the 
following conditions are satisfied: (i) the states fall into 
approximate multiplets of SU(2)l x SU(2)r and U(1)a 
and the splittings within the multiplets ( A/i = /i + — /i_) 
vanish at n — ► oo and/or J — > oo ; (ii) the splitting 
within the multiplet is much smaller than between the 
two subsequent multiplets. 

To verify this prediction we have to solve the Bethe- 
Salpeter equation, which will be done in the following 
sections. 



BETHE-SALPETER EQUATION FOR MESONS 
AND CANCELLATION OF THE INFRARED 
DIVERGENCES 

The homogeneous Bethe-Salpeter equation (BSE) for 
a quark-antiquark bound state with total and relative 
four momenta P and p, respectively, can generally be 
written in covariant form as 

X(P,P) = ~ij ^K(P,p,q)S(q + P/2) 



(2nY 

X X (P,q)S(q-P/2), 



(29) 



where K(P,p,q) is called the Bethe-Salpeter kernel and 
x(P, p) is the meson vertex function. In the rest frame, 
i.e. with the four momentum P v = (/x, P = 0) the latter 
depends on the meson mass /i and all four components of 
p. Finally, for the instantaneous interaction of our model 
it is energy independent and BSE becomes 



d 4 q 



(2tt)4 

x x{^,0)S{q - ///2,fe)7 . 



V{\p-<R) 10 S{q +n/2,k) 



(30) 



In order to solve the Bethe-Salpeter equation for 
mesons with the set of quantum numbers J PC with the 
given total spin J, parity P, and C-parity C, we expand 
the meson vertex function XjaA^iP) into a se t of all pos- 
sible independent Poincare-invariant amplitudes consis- 
tent with J PC . This is done in the Appendix A. Then 
the Bethe-Salpeter equation transforms into a system of 
coupled integral equations. Such a system of equations is 
derived in the Appendix B. The expansion of the meson 
vertex functions as well as the system of coupled equa- 
tions look slightly differently for three possible categories 
of qq mesons. 

Category 1. To this category mesons with J h for 
J = 2n and J H for J — 2n + 1 belong. In the non- 
relativistic limit these mesons have wave functions with 
orbital angular momentum L = J and the sum of the 
two quark spins S = 0. 

Category 2. Mesons with J++ for J = 2n and J~~ 
for J = 2n + 1 are in this category. In this case the 
nonrelativistic wave functions have L = J ± 1 and S = 1 . 

Category 3. This category contains mesons with J~ 
for J = 2(n + 1) and J++ for J = 2n+ 1. Their nonrel- 
ativistic wave functions have L — J and 5 = 1. 

The relativistic quark-antiquark states with J PC = 
do not exist in the BSE framework with an instan- 
taneous interaction. In the nonrelativistic limit L = 
and 5 = 1 would have to couple to J — 0, which is im- 
possible. The quark-antiquark states of category 4, i.e., 
with for J = 2n and J ^ for J = 2n + 1, do not 
exist in the BSE framework with the instantaneous in- 
teraction, too, as is also shown in Appendix B. There 
are also no corresponding nonrelativistic states. Mesons 
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with such quantum numbers are called exotic (mesons 
with exotic quantum numbers) and can be realized in 
form of hybrids, i.e., as qq states with explicit excitation 
of the gluonic field, or as four-quark states. If they can 
also be described as quark-antiquark states in the the 
BSE framework with a non-instantaneous interaction is 



a matter of dicsussion. We do not consider them in the 
present paper, however. 

Now we want to study the infrared limit of the Bethe- 
Salpctcr equation and cancellation of the infrared sin- 
gularities. To this end let us take a closer look at the 
integral equations for mesons of category 1, which are 



u)(p)h(p) 

g{p) = h{p) 



(31a) 



1 /' ,f'a M(p)M(q)Pj{p-q)+pq(0 I Pj +1 (p-q) + ^ I P J _ 1 {p-q)) 

w 9{qX31h) 



2 / (2tt) 3 



r 



Here the functions h(p),g(p) are as defined in the Ap- 
pendix B which along with the meson mass fx have to be 
found from the given equations and k = p — q. 

According to (|24[) and (121|) there are IR divergent 
terms 



2/ilR 



h(p) 



(32) 



and 



\ d 3 qSik)Pjip ■ q)hiq) = ^-hip) (33) 
2/im J 2^i R 



The term with giq) in the integrand of (|31ap comes with 
a factor l/w(<jf) = O(^ir) and the integral over this term 
becomes the IR finite expression 



-9ip)- 



Similarly there are IR divergent terms 



2mir 



lip) 



(34) 



(35) 



on the left and right hand sides of (|31ap . respectively. and 



Mip)Miq)Pjip -q)+pq (fi^Pj+lip ■ q) + ^Pj-xip • qj) 



CT 



., , d'qSik) Vrn = o 9{P) (36) 

2//ir J u}(p)0Jiq) 2/xir 

on the left and right hand sides of (|31b[) , respectively. The second term in the bracket on the left hand side has a 
factor cj(p) _1 = 0(/Xir) which vanishes in the IR limit. On the right hand side there is an IR finite term h(p). One 
sees that the IR divergent terms on both sides of both equations cancel and we end up with the coupled system of 
integral equations 



w/(p)fc(p) = ^9(P) + 1 J -^V f {k)Pjip-q)h(q), 

ufip)gip) = Hp) 



(37a) 



1 



d 3 q Mip)Miq)Pji P -q)+pq (^Pj+liP • q) + J^Pl-iiP ' <?)) 



Vfik) V"wl ■ ( 37b ) 



2 J (2tt)3 Jy ' w(p)w(g) 

with 



V f ik)=Vik) i2n) 3 Sik), (38) 

Mir 



which has a finite IR limit. 
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In an analogous manner all IR divergences from the integral equations for the mesons of category 2 and 3 can be 
removed. The IR finite integral equations for mesons in category 2 are 



^f(p) h i(p) 



9i(p) 



q 



V f (k) 



J 



2J+1 



Pj+tip-q) 



J 



2J+1 



M(q) y/J(J+l) 

Q(q) 2J+i 



{Pj+l{p-q)-Pj-l{p-q))h 2 {q) 



Pj-xip-q) h x {q) 



(39a) 



Uf(p)gi(p) 



d 3 q \ pqPj(P ■ q) + M (p)M (q) (jj^Pj+i (p • q) + ^Pj-liP • q)) 

Vf{k) \ ^WM) 



(27T) 



9%(q) 



M( P ) y/j(J+l) 

Q{p) 2 J + 1 



(Pj+x{p-q)-Pj-l(p-q))g2{q) 



(39b) 



u) f (p)h 2 (p) = 



— 92(P) + 2 



d 3 q (pqP.l(p ■ q) + M(p)M(q) fj±^P J+1 (p . $) + ^Pj^p ■ q)) 

V f( k ) \ ^Cw3 M<7) 



(27T) 



w(p)w(g) 



M(p) yj(j + i) 

w(p) 2J+1 



(Pj + l(p-g)-Pj-l(0-$))fcl(?) 



(39c) 



Vf(p)a2(jp) = h 2 (p) 



d 3 q 
(2^3 



/ J + l 
\2J+1 



Pj+i(p-q) 



+ 



M(q) y/J(J + l) 
w(g) 2J+1 



J 



2J+1 



Pj-xip ■ q) 92(q) 



(P J+1 (p-q)-Pj-l{p-q))9l{q) 



(39d) 



Note that the number of the auxiliary functions hi, <?i, h 2 , 92 as well as the number of coupled equations in the 
present case is twice larger than for mesons of category 1. This is related to the number of independent components 



in Xjm(p>P) f° r J > 0- However, for the states ++ within the same category only the last two equations (|39c 
(I39dp apply, because in this case there are only two independent functions h 2 , g 2 ■ 
For mesons in category 3 we have: 



Uf{p)h{p) 



-g(p) 



+ 



Uf(p)g(p) = h( P ) + - j 



d 3 q M(p)M(q)Pj( P -q)+pq (jj^P.r+lip ■ q) + ^jPj-iip ■ q)) 
3 V f( k ) ^T7Zw3 -%)' 

^V f (k)Pj(p-q)g(q). 



(27T) 
1 



j(p)uj{q) 



(40a) 
(40b) 



From the integral equations the normalization of the 
vertex functions cannot be determined. One can derive 
a normalization condition by demanding that the charge 
of the isovector mesons with isospin projection +1 must 
be one, or equivalently by normalizing the wave function 



by 



tr 



d 3 p 
(2tt) 3 



(^mHp,P)^ C jm{^P) 

-i, P - C JM\^i> P -°JMM)=^ (41) 



(see the Appendix C for the definition of the wave func- 
tion), where traces in Dirac space and color space have 
to be taken, the latter giving just a factor 3. Then for 



8 



mesons of categories 1 and 3 this yields 

I J Wf h ^P) = I- (42) 

For mesons of category 2 the normalization condition be- 
comes 

h^t j ^MP)*(P)=M+M = 1, (43) 

Here A/i and A/2 represent the normalization factors of 
two coupled amplitudes, (h±,gi) and (/i2,ff2)- Since 
these amplitudes are mixed by the dynamical mass of 
quarks, M(j>), a mixing angle olm can be defined by 

I 



M\ — sin 2 ttif and A/2 = cos 2 o^/ . Obviously in the limit 
M(p) — > the mixing vanishes. 



EFFECTIVE RESTORATION OF CHIRAL 
SYMMETRY FOR MESONS 

In this section we demonstrate that in the limit of the 
vanishing dynamical mass M(p) — > the Bethe-Salpeter 
equation of the previous section recovers all the antici- 
pated chiral multiplets (J25])- (J27J) . 

For M(p) = the integral equations (|37|) . which de- 
scribe mesons of category 1, i.e. , with J for J = 2n 
and J" 1 for J = 2n + 1, become 



M 2 _^ , I [ d 3 q 
2 J W) 



u f (p)h(p) = ^9(P) + u / 7^3 ' (44a) 



uf(p)g(p) = ft(p) + ^ y (2^)3 ^( fc ) 27TI ^ 

For mesons of category 2, i.e., with J ++ for J = 2n ( J 7^ 0) and J for J = 2n + 1, the four coupled integral 
equations (|39[) decouple into two independent systems of equations 



2 



1 f d 3 q T/ n JP J+1 (p-q) + (J+l)Pj_ 1 {p-q) 



wftohtw = ^ 51 ( P ) + - y ^-^(fc) "— ^ zj + 'i^"" 1 ^ " Mg)i (45a) 

L0 f (p) 9l (p) =h l (j ? )+ 1 -J ^Lv f (k)Pj(p ■ q) 9l {q) (45b) 



and 



W/(P)/» 2 (P) = ^U(p) + ^ / ^3 ^( fc )^(^ ' ( 45C ) 

n n , / x , 1 f d?q (J + l)P J+1 (p ■ q) + JPj-i(p ■ q) 
^f(j>)92{P) = h 2 [p) + - J J^ v f( k ) 2jTl ff2 ^' *• ' 

For the ++ states within this category only equations (|45c[) - (|45d|) apply. 

For the mesons of category 3, i.e., with J ++ for J = 2n + 1 and for J = 2(n + 1), the Bcthc-Salpeter equation 
takes the form 

r \ur \ V 2 f s , 1 [ d 3 q JP J+1 (p ■ q) + (J +l)Pj_i(p • q) 
Uf(P)h(p) = -j9(P) + 2 J j2nf Vf ( k > 2JT1 ^' ^ 



1 f d 3 



(1 



Uf(p)9(p) = Hp) + -J —?-V(k)Pj(p ■ q)g(q). (46b) 

I 



We see that the equations (|44a[) - (|44bj) become identi- cal to the equations (|45c|) - (I45dj) as well as the equations 
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(I45ap - (|45b[) are identical to equations (|46ap - (|46b|) . 
Note, that there is no isospin dependence of the interac- 
tion, i.e. each system of coupled equations describes both 
1=0 and 1=1 states which are exactly degenerate. Then 
we recover that the equations above fall into chiral mul- 
tiplets ([2S])-(I23 and both SU{2) L x SU{2) R and U{l) A 
are manifest. Hence if the typical momentum of quarks 
is large, as it is anticipated in the highly-excited mesons, 
the spectrum should be organized into chiral and U(\)a 
multiplets , because at large momenta dynamical mass 
of quarks M(p) vanishes. 



NUMERICAL RESULTS FOR SPECTRUM 

The method of solution for the BSE used in this work 
is described in Appendix D. In Table 1 we present our 
results for the spectrum for the two-flavor (u and d) 
mesons in the chiral limit m — 0. Due to chiral sym- 
metry breaking there is a mixing between (hi,g±) and 
(^2) 92) amplitudes for mesons of category 2 (but not for 
the ++ states). In the language of chiral representations 
this mixing means that for mesons of category 2 with 
J = 2n + 1 the representations (0, 0) and (1/2, l/2) a are 
mixed in the mesons 0, J as well as the representations 
(0,1) © (1,0) and (1/2,1/2);, are mixed in the mesons 
1, J~~. For mesons of category 2 with J = 2n, J > the 
representations (0,0) and (l/2,l/2) a are mixed in the 
0, J ++ states and the representations (0, 1) (1,0) and 
(1/2, 1/2);, are mixed in the mesons 1, J ++ . The mesons 
of the category 2 in the Table are assigned to a defi- 
nite chiral representation according to the chiral Bcthc- 
Salpeter amplitude which dominates in the meson wave 
function, which is determined by the relative magnitude 
of the normalization factors N\ and A2 in (|43|) . 

Note that within the present model there are no vac- 
uum fcrmion loops. Then since the interaction be- 
tween quarks is isospin-independent the states with the 
same J but different isospins from the distinct mul- 
tiplets (l/2,l/2) a and (l/2,l/2)b as well as the states 
with the same J PC but different isospins from (0, 0) 
and (0,1) © (1,0) representations are exactly degener- 
ate. Hence it is enough to show a complete set of the 
isovector (or isoscalar) states. 

The presented results are accurate within the quoted 
digits at least for states with small J and for states with 
higher J but small n. At larger J for larger n numerical 
errors accumulate in the second digit after comma. 

The axial anomaly is absent within the present model, 
because there are no vacuum fermion loops here. Even so 
there are no exact U (1)a multiplets, because this symme- 
try is broken also by the quark condensate of the vacuum. 
The mechanism of the U{\)a breaking and restoration is 
exactly the same as for SU(2) L x SU(2) R . 

The excited mesons fall into approximate chiral and 
U(1)a multiplets and all conditions of the effective sym- 



TABLE I: Masses of isovector mesons in units of ^/a. 



chiral 






radial excitation 


n 




multiplet 





1 


2 


3 


4 


5 


6 


(1/2,1/2). 


Q-+ 


0.00 


2.93 


4.35 


5.49 


6.46 


7.31 


8.09 


(1/2,1/2)6 


0++ 


1.49 


3.38 


4.72 


5.80 


6.74 


7.57 


8.33 


(1/2,1/2). 


1+- 


2.68 


4.03 


5.16 


6.14 


7.01 


7.80 


8.53 


(1/2,1/2)6 


1"" 


2.78 


4.18 


5.32 


6.30 


7.17 


7.96 


8.68 


(0,1)0(1,0) 


1~ 


1.55 


3.28 


4.56 


5.64 


6.57 


7.40 


8.16 


(0,1)0(1,0) 


1++ 


2.20 


3.73 


4.95 


5.98 


6.88 


7.69 


8.43 


(1/2,1/2). 


2" + 


3.89 


4.98 


5.94 6.80 


7.59 


8.31 


8.99 


(1/2,1/2)6 


2++ 


3.91 


5.02 


6.00 


6.88 


7.67 


8.40 


9.09 


(0,1)0(1,0) 


2++ 


3.60 


4.67 


5.64 


6.51 


7.32 


8.06 


8.75 


(0,1)0(1,0) 


2" 


3.67 4.80 


5.80 


6.68 


7.49 


8.23 


8.91 


(1/2,1/2). 


3+" 


4.82 


5.77 6.62 


7.41 


8.13 


8.81 


9.45 


(1/2,1/2)6 


3"" 


4.82 


5.78 


6.65 


7.44 


8.17 


8.86 


9.50 


(0,1)0(1,0) 


3~ 


4.68 


5.63 


6.48 


7.26 


7.99 


8.67 


9.30 


(0,1)0(1,0) 


3++ 


4.69 


5.66 


6.53 


7.32 


8.06 


8.74 


9.39 



(1/2,1/2). 


4~ + 


5.59 


6.45 


7.23 


7.96 


8.64 


9.28 


9.89 


(1/2,1/2)6 


4++ 


5.59 


6.45 


7.24 7.97 


8.66 


9.30 


9.91 


(0,1)0(1,0) 


4++ 


5.51 


6.36 


7.15 


7.88 


8.56 


9.19 


9.80 


(0,1)0(1,0) 


4— 


5.51 


6.37 


7.16 


7.90 


8.58 


9.23 


9.84 


(1/2,1/2). 


5+" 


6.27 


7.05 


7.78 


8.47 


9.11 


9.72 


10.3 


(1/2,1/2)6 


5"" 


6.27 


7.06 


7.79 


8.47 


9.12 


9.73 


10.3 


(0,1)0(1,0) 


5~ 


6.21 


7.00 


7.73 


8.41 


9.06 


9.67 


10.2 


(0,1)0(1,0) 


5++ 


6.21 


7.00 


7.73 


8.42 


9.07 


9.68 


10.3 


(1/2,1/2). 


6-+ 


6.88 


7.61 


8.29 


8.94 


9.55 


10.1 


10.7 


(1/2,1/2)6 


6++ 


6.88 


7.61 


8.29 


8.94 


9.56 


10.1 


10.7 


(0,1)0(1,0) 


6++ 


6.83 


7.57 8.25 


8.90 


9.51 


10.1 


10.7 


(0,1)0(1,0) 


6~ 


6.83 


7.57 


8.26 


8.90 


9.52 


10.1 


10.7 
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FIG. 3: Spectra of J = mesons (masses in units of 



metry restorations are satisfied. We observe a very fast 
restoration of both SU (2)l x SU(2) r and U(1)a symme- 
tries with increasing J and essentially more slow restora- 
tion with increasing of n. The spectrum of the J = 0, 1, 2 
mesons is shown in Figs. [3] -El 

In Fig. [S]the rates of the symmetry restoration against 
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FIG. 6: Mass splittings in units of y/a for isovector mesons 

S ~ 3 of the chiral multiplets (1/2, l/2) and (1/2,1/2)6 (circles) 

and within the multiplet (0, 1) © (1, 0) (squares) against J for 

|__ | n = (top) and against n for J = and J = 1, respectively 

o,i+- 1,1" 1.1++ (bottom). The full line in the bottom plot is 0.7^/a/n. 



FIG. 4: Spectra of J — 1 mesons (masses in units of i/a). 



(0:0) (1 
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(0,1)0(1,0) 



- L0 

9 

- 8 



L2++ 



FIG. 5: Spectra of J = 2 mesons (masses in units of 



the radial quantum number n and spin J are shown. It is 
seen that with the fixed J the splitting within the multi- 
plets AM decreases asymptotically as 1/ \/n, dictated by 
the asymptotic linearity of the radial Regge trajectories 
with different intercepts. With the fixed n the J-rate of 
the symmetry restoration is much faster. 

In the limit n — > oo or J — > oo one observes a complete 
degeneracy of all multiplets, which means that the states 
fall into representation (|28|) that combines all possible 
chiral representations for the systems of two massless 
quarks [7(. This means that in this limit the loop effects 
disappear completely and the system becomes classical 

d, Q. The analytical proof for this larger degeneracy 
will be given in the next section. 

In Fig. [7] the angular and radial Regge trajectories 
are shown. Both kinds of trajectories exhibit asymptot- 
ically linear behavior. This has to be expected a-priory, 
because at large J or n all higher Fock components are 
suppressed and asymptotically vanish (it will be later well 
seen from the meson wave functions). Then the semiclas- 
sical description of the states with large n or J with the 
linear potential requires the linear Regge trajectories, see, 

e. g., ref. [3§|. There are deviations from the linear be- 
havior at finite J and n, however. This fact is obviously 
related to the chiral symmetry breaking effects for lower 
mesons. Note, that the chiral symmetry requires a dou- 
bling of some of the radial and angular Regge trajectories 
for J = 1,2, ... This is a highly nontrivial prediction of 
chiral symmetry. For example, some of the rho-mesons 
lie on the trajectory that is characterized by the chiral 
index (0,1) + (1,0), while the other fit the trajectory with 
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n 



FIG. 7: Angular (top) and radial (bottom) Regge trajectories 
for isovector mesons with M 2 in units of a. Mesons of the chi- 
ral multiplet (1/2, l/2) a are indicated by circles, of (1/2, l/2)& 
by triangles, and of (0, 1) © (1,0) by squares (J ++ and J~~ 
for even and odd J, respectively) and diamonds (J and 
J ++ for even and odd J, respectively). 



the chiral index (1/2, l/2) 6 . 

The numerical result for the quark condensate is (qq) = 
(— 0.231v / o : ) 3 , which agrees with the previous studies 
within the same model. If we fix the string tension from 
the phenomenological angular Regge trajectories, then 
\fo w 300 - 400 MeV and hence the quark condensate is 
between (-70 MeV) 3 and (-90 MeV) 3 which obviously 
underestimates the phenomenological value. Probably 
this indicates that other gluonic interactions could also 
contribute to chiral symmetry breaking. 



Let us compare the integral equations (|37| and (|40f for 
mesons of categories 1 and 3, respectively. By exchanging 
h and for the latter (|40b[) becomes (|37a[) and the 
difference between (|37b[) and (|40ap reduces to different 
coefficients in front of the Legendre polynomials, i.e. one 
has 

(J+l)Pj+i(p-g) + JPj-i(p-q) , 47 , 
2J+1 1 ' 

and 

27TT ' (48) 

respectively. The same factors are also seen in eqs. (|44a|l - 
(|46bp . The difference between (gTJ) and (gSJ) vanishes 
for large J. This means that for large J there appears 
degeneracy within the whole reducible chiral multiplet 
[(0,1/2)0(1/2,0)] x [(0,1/2)0(1/2,0)] which combines 
all possible chiral representations of the quark-antiquark 
systems with the same J. 

Will we see similar higher degree of degeneracy at large 
radial quantum number n but finite J? The answer is yes, 
which is demonstrated in what follows. 

In general, for a large meson mass, i.e. n and/or J 
large, the quark-antiquark component ip + of the meson 
wave function propagating forward in time be dominant 
against the backward-propagating component ^_ (see 
Appendix C). In the IR limit n/w(p) in (fW|) - (fTTJU]) goes to 
zero, which means that for all three categories of mesons 
4> + ~ l/2(h + %g) and ip_ ~ 1/2 (ft - %g). For highly 
excited mesons, by neglecting the latter and setting 

'0+ » h « fig, (49) 

one ends up with only one equation for ip + for all quark- 
antiquark states with the given J and n, i.e., 

w/(p)^+(p) = 2^+(P) 

Hence all possible states with the same J at the given 
large n fall into the [(0,1/2) © (1/2,0)] x [(0,1/2) © 
(1/2,0)] representation and the degeneracy becomes ex- 
act when J and/or n approach infinity. 



HIGHER REPRESENTATIONS FOR EXCITED 
STATES 

We have mentioned in the previous section that an 
approximate degeneracy within the higher representation 
(|28p is observed for highly excited mesons. Below we 
demonstrate this property analytically. 



WAVE FUNCTIONS 

At large J or n the semiclassical description requires 
that the higher quark Fock components be suppressed 
relative the leading one and asymptotically vanish. This 
is because the higher Fock components manifestly repre- 
sent effects of quantum fluctuations. The leading Fock 
component is the forward-propagating quark-antiquark 
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FIG. 8: Components ip+(p) and ip-(p) of the wave functions 
for ground and radially excited states for mesons of cat- 
egory 1. All quantities are given in appropriate units of y/a. 
Notice that the ground state (Goldstone boson) has /i = 
and thus V'-(p) = V*+(p)j hence the two curves for n = 
coincide. 
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FIG. 9: Components ip+(p) and ip-{p) of the wave functions 
for mesons J~ + of category 1 with radial quantum number 
n = and spin J — 2, 4, 6. All quantities are given in appro- 
priate units of y/a. 

component tp+(p), defined in the Appendix C, while 
the higher Fock components contain necessarily the 
backward-propagating quark- antiquark states tp-{p)- I n 
Figs. [8] and [9] we show that for mesons of category 1 
with large radial quantum number n and/or spin J the 
forward-propagating component ijj+(p) becomes indeed 
much larger than the backward-propagating component 
^_(p). The wave functions of mesons of categories 2 
and 3 behave similarly. In order to study the behavior of 
the wave functions of mesons with large n and/or J we 
can thus restrict ourselves to ip+ (p) . 

In Fig. rrU] we compare the wave functions of chiral 



FIG. 10: Components of ground states and radial exci- 

tations of pseudoscalar and scalar mesons. All quantities are 
given in appropriate units of y/a. 



partners in (1/2,1/2) representations for mesons with 
^ and ++ , for the ground states (n = 0) and radial 
excitations n = 3, 6. A principal difference between pseu- 
doscalar and scalar mesons is the behavior of the wave 
functions for momentum p — > which go to a constant for 
the former and to zero like p for the latter. Thus at small 
p there is a large difference between the wave functions 
of pseudoscalar and scalar mesons with the same radial 
quantum number n even if n is large. At large p (where 
dynamical mass becomes small and effective restoration 
of chiral symmetry takes place) the wave functions of the 
radially excited pseudoscalar and scalar mesons behave 
similar but there is a shift in the positions of nodes and 
maxima/minima. This persisting difference of the wave 
functions is a reason for a slow rate of chiral restoration 
in meson masses with increasing n. 

For J > there are four different mesons with the 
given isospin at fixed n and J. One of category 1 and 3, 
respectively, and two of category 2. 

Consider mesons of category 2. They are mixings of 
two different chiral representations with radial functions 
■01+ (p) and il , 2+(p)i respectively The mixing is charac- 
terized by the size of the normalization factors J\f\ and 
Af 2 = 1 In Fig. [TT] we show A/i (A/2) for all states 

of category 2 considered in this work. For given n the 
mixing decreases with increasing J. At large J there is 
one almost unmixed state with A/i. ~ 1 and one almost 
unmixed state with A/2 w 1. The former belongs to the 
chiral multiplet (0, 0) and (0, 1) (1, 0) in case of isospin 
and 1, respectively and the latter to (4, |) a and (|, |)b 
in case of isospin and 1, respectively. In contrast, for 
given J the mixing first increases with increasing n and 
then saturates. The saturation value seems to be of about 
0.6 for J = 1. For larger J such a saturation is not yet 
reached for n — 6. 
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FIG. 11: A/i (or A/2 = 1 — Afi) for mesons of category 2. For 
each J > and n there is one state with Afi > A/2 = 1 — A/i 
and vice versa. 
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FIG. 12: Components ip+(p) of mesons with J = 1. Upper 
plot: n — 0, lower plot: n = 6. All quantities are given in 
appropriate units of 



In Fig. [12] we show the wave functions For J — 1 and 
n = and n = 6. In both cases the wave functions 
of mesons of category 1 and 3 are rather similar to each 
other. These wave functions go to zero for p — > contrary 
to the wave functions of the mesons of category 2 which 
go to a constant. At small momenta the coupling between 
ipiip) and ip2(p) m mesons of category 2 becomes large 
and thus there is a rather strong mixing even for high 
radial excitations. 

The radial wave functions of states with n = and in- 
creasing J are shown in Fig. [T3]for J = 2,4,6. Already 
for J — 2 four of these functions look rather similar. 
These are the ones for the mesons of category 1 and 3 
and in each case the larger function of the two mesons of 
category 2, which is in the first case and in the 

second case ip2+(p)- The two other functions (?/>2+(p) 
and in the first and second case, respectively) are 

smaller. From Fig. [Til one sees that Afi (or A/2) are both 
larger than 0.8 in the first and second case, respectively. 
For J = 4 and even more for J = 6 the mixing becomes 
very small and consequently the physics is determined by 
the four almost indistinguishable radial functions. Con- 
sequently there appears approximate degeneracy of all 
mesons for given n at large J. This is in accord with the 
discussion in the preceeding section. Obviously the wave 
function for larger J are strongly suppressed at small 
momenta, i.e., the radial function goes to zero as p J for 
mesons of categories 1 and 3 and as p ,7_1 for mesons of 
category 2. It is this feature which explains a fast chiral 
restoration with increasing J, because the chiral symme- 
try breaking dynamical mass M(p) is essential only at 
small momenta. 

For higher radial excitations the wave functions ac- 
quire nodes and the first bump is shifted towards lower 
momenta as we show in Fig. [14] for J — 6 and n = 6. 
That explains why at given J for increasing n the mixing 
becomes larger. 



CONCLUSIONS 

In this paper we have demonstrated explicitly a fast 
effective chiral symmetry restoration in excited mesons 
with increasing J and a slow restoration with increasing 
n. To this end we used a model where the only interac- 
tion between quarks is linear instantaneous Coulomb-like 
confining interaction between color charges. This model 
contains all the necessary ingredients - it is confining and 
manifestly chirally symmetric. Chiral symmetry break- 
ing is provided in the standard way through the non- 
perturbative self interactions of quarks while the meson 
spectrum is obtained from the Bethe-Salpeter equation. 

The key implication of dynamical chiral symme- 
try breaking is that the quarks acquire the dynami- 
cal Lorentz-scalar mass which is exclusively of quan- 
tum (loop) origin. This dynamical mass is strongly 
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FIG. 13: Components ip+(p) of mesons with J = 2, 4, 6 and 
7i = 0. All quantities are given in appropriate units of y/a. 



momentum-dependent and vanishes fast at larger mo- 
menta. In highly excited states the quantum (loop) 
effects must be suppressed relative the classical con- 
tributions and the spectrum exhibits an approximate 
SU(2)l x SU(2)r and U(\)a symmetry. Microscopically 
this is because a typical momentum of valence quarks in- 




FIG. 14: Components ip+(p) of mesons with J = 6 and n = 6. 
All quantities are given in appropriate units of y/a. 



creases higher in the spectrum and consequently their dy- 
namical Lorentz-scalar mass, which violates chiral sym- 
metry, becomes small and asymptotically vanishes. 

Within nonrelativistic or semirelativistic quark mod- 
els the effective chiral restoration cannot be reproduced 
as a matter of principle, because (i) the constituent 
mass of quarks is an input parameter and is momentum- 
independent, (ii) the notion of the quark chirality is ab- 
sent within a nonrelativistic or semirelativistic approach, 
(iii) a Lorentz-scalar confining interaction, which is typi- 
cally assumed in these quark models, breaks chiral sym- 
metry manifestly. 

We have demonstrated explicitly that in the high- 
lying mesons all higher Fock components are strongly 
suppressed relative the leading qq one and a very fast 
restoration of chiral symmetry happens with increasing 
J and rather slow one with increasing n. A reason for 
the former is that with larger J the radial wave function 
vanishes at small relative momenta, and hence the large 
chiral-symmetry breaking dynamical mass becomes irrel- 
evant. The slower rate of the restoration in the latter 
case is due to the fact that at J = the wave function at 
small momenta does not vanish for pseudoscalar mesons, 
so that even though a typical momentum of quarks in- 
creases with n, the small momentum dynamical mass of 
quarks always contributes to some extent. 

All possible states for a given n and large J fall into 
higher representations [(0,1/2) © (1/2,0)] x [(0,1/2) © 
(1/2,0)] that combine all possible chiral multiplets with 
the given J and n. The same property is observed for 
large n at a given J, though a degree of degeneracy is 
smaller. This property means that in both these cases 
the quantum loop effects become irrelevant and all pos- 
sible states with different orientations of quark chirality 
become equivalent. 

We have shown that both radial and angular Regge tra- 
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jectories are linear asymptotically. At smaller J and/or n 
the linear form is broken due to chiral symmetry break- 
ing effects. Each Regge trajectory is characterized by 
the proper chiral index and hence the amount of inde- 
pendent Regge trajectories coincides with the amount of 
chiral representations. Consequently for many mesons 
there are two independent Regge trajectories at given 
/, J PC , because in the chirally restored regime there are 
two independent and degenerate states with the same 
/, J which belong to different chiral representations. 



ACKNOWLEDGMENTS 



RFW acknowledges helpful discussions with R. 
Alkofer, M. Kloker, and A. Krassnigg. This work 
was supported by the Austrian Science Fund (projects 
P16310-N08 and P19168-N16) and by DFG (project Al 
279/5-1). 



J = 2n + 1, respectively have the vertex function 

+ i lfl e^(el J ' M) (P) ■p)P p p a {p ■ p, 

+ i^pM jm \p)-p)xz 

+ 1sYpM JM) (P) ■ P)(P ■ P)X4 

+ i 5 r(^' M \p) ■ p)(p ■ P)x 5 

+ ii 5 ^ v p^M jm) (p)-p)xs 

+ ll^ V p^i J - M \P) ■ P )(p • P) XS - 



Here 



X.2 



(e« M )(P).p) = e« M j J (Pr.../ J , 



(51) 



(52) 



(#' M) ( p ) ■ P) = i P )P" 2 ■ ■ ■ P" J > (53) 



and 



MM) 



j (P) are the polarization tensors of rank J. 



They have the properties 



MM) 



(54) 



9 fHH ^fL,...(P) = 0, 



APPENDIX A 

In this Appendix we discuss a general parametrization 
of the vertex function and its specific form for mesons of 
all categories with the instantaneous interaction. 

For a meson with mass /i, spin J, intrinsic parity P 
and C-parity C the vertex function Xjm(P-iP) can be 
written as a sum of eight (four for J = 0) components, 
where M = —J,..., J is the z-projection of J in the 
rest frame of the meson. P v and p v are the total and 
relative four-momenta of the quark-antiquark system, re- 
spectively. Each component contains a function x% which 
depends on the two invariants p 2 and (p.P) 2 (notice that 
p.P is Lorentz-invariant but changes its sign under the 
C-parity transformation) . Since we are only interested in 
spectra in this paper we can restrict ourselves to the rest 
frame of the meson (i.e., the frame with P v = (/x, 0, 0, 0)), 
where the vertex function depends on fi, pq and p and 
each component contains a function Xi °f the two vari- 
ables Pq and \p\ . Some of the components contain a factor 
Po and are thus odd functions of po- At the same time, 
for an instantaneous interaction the vertex function in 
the rest frame is independent of po. Thus all components 
which are odd in po must vanish. The function Xi °f the 
remaining components depend only on p = \p\. Depend- 
ing on the quantum numbers J PC there are four different 
categories of mesons. In the following we give the general 
expressions for the vertex function X^MiP^P) m a general 
frame, XjAtil^'POiP) m the rest frame and XjAii^'P) m 
the rest frame and an instantaneous interaction. 

Category 1: Mesons with J ^ for J — 2n and for 



and 



e~r >{P)P f2 



MM) I 

.41.. 



0. 



(55) 



(56) 



Contracting a polarization tensor with J vectors at, i = 
1, ... J yields in the rest frame 

e K ) J ((/x ,0,0,0))af ■■■a^ = 

j (2J—i)\\ \ ) 

y K j\ — {{■ • . {ai <g> a 2 } 2 ® ■ ■ -}j-i O a,j}j M , 

where { J\ ® J-i\jm is the usual coupling of two spherical 
tensors of rank J\ and J2 to a spherical tensor of rank J . 
Thus the vertex function (1511) can be written in the form 



Xjm(m,-Po,p) = 75*jm(p)xi(Po,|p1) 

+ M7o7s50m(p)X2(Po! IpI) 
+ M7075 Oj+ 1 (p) ® 7} ,7 M X3 (Po , \P\ ) 
+ M7o75 {Yj-i (p) O 7} j M X4 (Po , \p\ ) 
+ Pq {Yj(p) ® j}jm X5 
+ PoM7o {Yj(p) ® 7} JM XeiPo, \p\) 
+ Po^15{ y J+i{P)®i}jmX7{pIAp\) 
+ PQV>75{Yj-i(p)®j} JM Xs(Pa,\p\), 

(58) 

where the Yjm{p) are the spherical harmonics. For J = 
the components 4, 5, 6, 8 are absent. The components 1- 
4 are even and the components 5-8 odd functions of po, 
respectively. For an instantaneous interaction we are left 
with the vertex function 

XjM( m ,P) = 75 Y JM (p)xi (p) + V1qJ5Yjm(p)X2(p) 



+ M7o75{y>+i(p) ® i}jmX3(p) 
+ M7o75 {X/-i(p) ® l} JM Xi{p) 



(59) 
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for mesons of this category. 

Category 2: Mesons with J ++ for J = 2n and for 
J = 2n + 1 , respectively have the vertex function 

x5&(p,p) - {t {JM) (P)-p)xi 

+ h 5 ^e^(ei J ' M \P)-p)P pPa X2 



(60) 



+ rP^ (J ' M \P)-p)(p-P)X3 

+ i»pM hM) {P)-p)xi 

+ r(ei JM) (P)-p)x5 

+ ia^P pPu (e^ M \P)-p) X 6 

+ ia^P^ J ' M \P)-p) X 7 

+ ia^Pp^ LM) (P)-p)(p-P) X8 . 



In the rest frame it can be written as 

Xjm(^Po,P) = Yjm{p)xi(pIAp\) 

+ {Yj+i(p)®1}jmX2{pI\p\) 

+ {Yj-i{p)®i}j M Xz{pI,\p\) 

+ Vl5{Yj{p)®^} JM Xi{plAp\) 

+ mo{Yj+i{p)®l} JM X5{pl,\p\) 

+ M7o {Yj-i{p) ® 7>j M Xefao, \p\) 

+ PoJoYjm(p)X7(pI,\p\) 

+ PoM7o75{^j(p)0 7}jmX8(Po-Ip1)- 

(61) 

For J = the components 3, 4, 6, 8 are absent. The com- 
ponents 1-6 are even and the components 7,8 odd func- 
tions of po, respectively. For an instantaneous interaction 
we are left with the vertex function 



Xjm(^P) = Yjm(p)xi(p) + {Yj + i(p)®i} jm X2{p) 

+ {y>-i(p)®^}jiifX3(p) 

+ M75 {Yj(p) ® 7}jm Xi{p) 
+ M7o {Yj+i{p) ® 7}jm X5{p) 
+ Wo{Yj-i(p)®i} JM X6(p), 

(62) 

for mesons of this category. 

Category 3: Mesons with J for J = 2n and J ++ for 
J = 2n + 1 , respectively have the vertex function 

x p j^{p,p) = i^ {j ' m) (p)-p){p-p)xi 

+ il^P°(ei JM) (P)-p)P p p a x 2 

+ i^pM JM) (p)-p){p-p)xz 

+ l5l^(e {J ' M) (P)-p)X4 

+ i^{^' m \p)-p)xb 

+ ii 5 <?^P,.pM JM) (P) ■ P)(P ■ P)X6 



In the rest frame it can be written as 

Xjm(m,Po,p) = v{Yj(p)®j} JM Xi(Po,\p\) 

+ l5{Yj+i{p)®l}j M X2{plAp\) 

+ 75{50-i(p)®7}jmX3(Po.|p1) 
+ lo{Yj(p)®7}j M X4(pl\p\) 
+ Po^YjmX5{pIAp\) 

+ POlOl5YjMX6(Po,\P\) 

+ Polol5 {Yj+-l(p) ®1}jmX7{pI,\P\) 

+ Po7o75 {Yj-i(p) ® 7>jm X&(Po, \p\) 

(64) 

For J = the components 1, 3, 4, 8 are absent. The com- 
ponents 1-4 are even and the components 5-8 odd func- 
tions of po, respectively. For an instantaneous interaction 
we are left with the vertex function 

Xjm(p,P) = v{Yj{p)®l} JM X\{p) 

+ 75{^/+i(p)0 7}j M X2(p) 
+ l5{Yj-i{p) ®l} JM Xz{p) 
+ 7o {Yj{p) <8> 7} JM Xi{p)- 



(65) 



for mesons of this category. 

Category 4: Mesons with J H for J = 2n and J h for 
J = 2n + 1, respectively have the vertex function 

X^&M - {e (JM) {P)-p){p-P)Xi 

+ «757^ AWpCT (^ J ' M) (- p ) ' P)i>«r(P • P)X2 

+ rpM ,M) (p)-p)x 3 
+ rp^ {J ' M) (p)-p)(p-p)x4 

+ 7n4 J ' M) (n-p)(p-^)X5 
+ i<J^P»Pv{e {JM) {P)-p){p-P)xs 
+ ta^P^ J ' M \p)-p)(p-P)x7 
+ in^p^ M \p)-p)x & . 

(66) 

In the rest frame it can be written as 

Xjm(v,Po,P) = M7o^jm(p)Xi(Po>Ip1) 

+ 7o75{ij(p)0 7}j M X2(Po>|Pl) 

+ po^y jm {p)xs{pIAp\) 

+ PoM^j+i(p)®7}jmX4(Po>Ip1) 
+ pom{^/-i(p) 7}j M X5(PoJp1) 
+ Po75{^/(p)0 7}j M X6(PoJp1) 
+ Po7o{^/+i(p)®7} jm X7(PoJp1) 
+ Polo{Yj-i(p) ®i} jm Xs,{pI, IpI) 

(67) 

For J = the components 2, 5, 6, 8 are absent. The com- 
ponents 1,2 are even and the components 3-8 odd func- 
tions of po, respectively. For an instantaneous interaction 
we are left with the vertex function 



+ il^P p {ei J ' M) {P)-p){p-P)x 7 
+ il^p^' M \P)-p)x 8 - 



Xjm(^,P) = H1oYjm{p)x\{p) 

+ 7o75 {Yj(p) ® j} JM X2(p) 



(68) 



(63) for mesons of this category. 
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APPENDIX B 

Given the parametrizations of the vertex function in 
Appendix A, we derive in the present Appendix from 
the general Bethe-Salpeter equation systems of coupled 
equations for all possible mesons categories. 



given quantum numbers becomes an equation between 
two Dirac-matrices. By projecting out the functions 
Xi(p) on the left hand side one arrives at a system of 
coupled integral equations. The projection can be per- 
formed by applying the standard rules of traces of Dirac 
matrices and with help of 



With the respective ansatz for the vertex function as 
given in the previous section the BSE for a meson of 



47T 



- J —Yj(p).Yj(q) = P J (p.q), 



(69) 



4tt 



2 J 



1 £{*Ji(p) ® e k }j ■ {Yj 2 {q) ® e k }j = (-ly+^Sj^Pj.ip ■ q), 



fe=i 



(70) 



47T 



2J + 



jYj(p)-{Yj 2 (q)®q}„ 



2J+l Oj201 ° 



Pj(p ■<?)=< 



J 



2J+1 
J+l 
2J+ 1 



Pj(p-q), J 2 = J-1 

Pj(p-q), J 2 = J+1 
0, else, 



(71) 



47T 



J 



2J+ 1 



Pj-i(p-q), Ji=J-l 



I±±P J+1 (p.q), J 1 = J+l 
0, else, 



(72) 



47T 



J+J VljTr^ o ^ 0oi0jPji (?3 
• 7 pj_i(p-?), j 1 = j 2 = j-i 



2J+ 1 



^^ Mtfg, J 1= J-1,J 2 = J+1 

^2J+ + l 1) pj+l( P ' ^ Jl = J + lj J2 = J ~ 1 
J+l „ , 



+ 
2J+ 1 



Pj+i(p-?), Ji = J 2 = J+l 
0, else, 



(73) 



47T 3 

j,k,l=i 



(-l) Jl+j2+ V6(2J 2 + l)C Jl ° 



/..in.i ^ ^ j f PjAv-q) 



= < 



J + i 

2J+ 1 



Pj-iip-q), Ji = J-1,J 2 = J 



^±lpj(p.g), J 1 = J,J 2 = J-1 



J 



2J+ 1 



fj(p-g), Ji = J,J 2 = J + i 



J 

2J+ 1 



Pj+i(p-g), Ji = J + l, J 2 = J 
0, else. 



(74) 
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For mesons of category 1 the BSE becomes a coupled system of four integral equations 
I 



1 f d 3 q V{k) 
XliP) = 2 J Wf^{^ Pj{lP - q) 



2Q(q) 



M(q)x2{q) - q 



J + l 
2J+ 1 



X3(q) 



J 



2J+1 



Xi(q) 



(75a) 



( s I f A V(k) 
X2{P > 2 J (2tt)3 w a (g) _ £ 



M( ? ) f 1 



4 



Xl(«) ■ -/ V 

2 



M(q)x2(q) - 9 



J+l 
2J+ 1 



X3(«) 



J 



2J+ 1 



X4(?) 



(75b) 



X3(P) = - 



J+l 1 f d A q V{k) 



2J+ 1 2 



(27r) 3 o; 2 (g) 



■Pj+i(p-q) 



w(q) 



1 / \ 



M(q)x2(q) ~ q 



J+l 
2J+1 



X3(q) 



J 



2J+1 



Xi(q) 



(75c) 



Xa(p) 



./ 1 f d 3 q V(k) 



P.J-i(p-q) 



2J+12J (2tt) 3 W 2(g) _ ii 

M(q)x2(q) - q 



q J 1 , , , w(g) 

0X1(9) 



w(g) 2' 



w(g) 



J + l 
2J+ 1 



X3(<?) 



J 



2J+ 1 



X4(g) 



(75d) 



Obviously in the integrands there appear only two linear combinations xi (q) & n d the term between the square brackets. 
Thus one can reduce (|75p to a coupled system of two integral equations 



uj(p)h(p) 



(76a) 



" £{pj) 9{P) = h{p) 



1 



2 J (2tt) 3 



d 3 g M(p)M(q)Pj( P -q)+pq (j±LP J+1 ( P ■ q) + ^Pj_ x (p . q)\ 

V(k) w r ^g(q](76b) 



w(p)w(g) 



for the two functions 



h(p) = 



Xi(p) 
uj(p) 



and 



g(p) 



h(p) + 2——X2{p) - 2 — 

uj(p) uj(p) \ V 2J + 1 



\/Stt X3(p) ~^jtt X4(p) 



(77) 



(78) 



r 



For J = Eq. (|75dp for ^4 is absent and all terms For mesons of category 2 the BSE becomes a coupled 
containing Xi m the integrands of (I75all75c| vanish. system of six integral equations 
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1 f d?q V(k) ,„ s 

Xlip) = 2 J wr^r^ Pj{p - q) 



Hi 

4 



qXi(q)+M(q) 



inh X2{q) V zjTi X3iq) 



+ ^ 
2 



27TT X5(9) -V2JTT X6(9) 



(79a) 



X2(P) 



i /■ d 3 9 y(fc) 

W ot , 1 V o T , 1 X2(9) 



+ 



+ 



J 

2J+ 1 

uj{q)M{q) 



J+l 



/i 



2 



J+l 



2J+ l^^' V 2J+ 1 
q X i(q)+M(q) I 



X3(q) 



J 



+ 



2Q(q) V 2J+ 1 
/x 2 M(g) 



2J+ 1 

g X 4(<z)+M(g) 



J+l 



2J+ 1 



X2(g) - 



J 



2J+ 1 



X3(g) 



J 



2J+ 1 



Xs(g) 



J+l 
2J+ 1 



Xe(g) 



J + l 



2u%) V 2J+ 1 



J+l 
2J+ 1 



xs(g) 



J 



2J+ 1 



Xe(g) 



(79b) 



X3(P) = 



d 3 q V(k) 



(27r) 3 w 2 (g) 



Pj-iip-q) 



J+l 
2J+ 1 

cj{q)M(q) 



[J: 



tf(q) 

2 



2J+ 1 

9Xi (?) 



X2(<?) + 



J+l 

2J+ 1 



xs(g) 



J + i 



2u(q) V 2J+ 1 
M 2 M(g) 



2J+ 1 

9X4(g) 



M(q) 
M(q) 



J 



2J+ 1 



J+l 
2J+ 1 

Xs(g) 



X2(g) 



2J + 



yXs(9) 



J+l 
2J+ 1 



Xe(q) 



J 



2Q{q) V 2J+ 1 



J+l f V 



J 



2J+ 1 



Xe(g) 



(79c) 



X4(P) 



d 3 g F(fc) 
(2tt)3 ^) _ 



-Pj(p ■ q) 



+ 



2w(g) 

qu(q) 

tf{q) 



j 



2J+ 1 



X2(g) + 



J+l 
2J+ 1 



X3(g) 



9X4(9) + M(q) (J^XM + /^Xeto) 



(79d) 
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X5(P) = 



1 r d 3 q V(k) 

2 J (2TTf L0 2 {q) _E 

fM(q) 



-Pj+lip-q) 



J 



[2uj(q)\ 2J + 1 [V 2J+1 

1 / J+l 
+ 2uj(q)\l 2 J + 1 

, M{q)uj{q) 



X2O2O 

q X l(q)+M{q) 



uj 2 (q) 
-ui(q) 



J 



2J+1 



J+l 
2J+ 1 



X3(g) 



J + i 

2J+1 



X2O2O - 



J 



2J+ l 



X3(g) 



qXi(q)+M(q) 



J 



2J+ 1 



xs(g) 



J + i 



2J+ 1 



J+l 
2J+ 1 



Xs(<?) - 



J 



2J+ 1 



Xe(g) 



J+l 
2J+ 1 



Xe(g) 



(79e) 



1 /■ d 3 g V(fc) 

X6(p) = 5/ (2^^r^ p ^-^ 



/ M(g) fl+T 



[2u;(g) V 2J+ 1 [V 2J+ 1 
1 



J 



J 



2w(g) V 2 J + 1 



X2(«) 

gXi(g)+M(g) 



J + l 
2J+ 1 



X3(q) 



J + i 

2J+ 1 



X2(g) 



J 



2J+ 1 



Xs(g) 



2J+ 1 



M(g)o;( g ) / J + l 
u) 2 (g) 

-w(g) 



gX 4 (g)+Af(g) 



J 



2J+ 1 



Xs(g) 



J 



2J+ 1 



Xe(g) 



J 



2J+ 1 



J+l 
2J+ 1 



Xs(<?) - 



J 



2J+ 1 



Xe(?) 



(79f) 



Obviously in the integrands there appear only four different linear combinations of the six Xi(p) given by the terms 
in the square brackets. Thus one can reduce (|T9[) to a coupled system of four integral equations 



w(p)hi(p) 



1 



V(k) 



d 3 q 

M(q) y/J(J + 1) 
ui{q) 2J+1 



P J+1 {p-q) + —^-Pj- 1 {p-q)) ih 1 {q) + -^gx{q) 



2J+ 1 



2J+ 1 



4w(q)' 



(Pj+l(p-q)-Pj-l(fi-q)) h 2 (q) 



4uj(q) 



52(g) 



(80a) 



w[p) 



4w(p) 



3i (p) = hi{p) 



cPq I PqPj(P ■ <l) + M(p)M(q) (^P J+1 (p ■ q) + Jd^Pj^p . 

3 V ( k ) S ~^7Iw3 L 9i{q) 



(27T) 



w(p)w(g) 



(80b) 



w(p)/ia(p) = 7: 



(27T) 

x ( /i 2 (g) 



d 3 ? j pgPj(0 • g) + M(p)M(q) (j±Lp J+1 (p ■ g) + -J—Pj_^ . q) 



u(p)u(q) 



M?) 



32(g) 



Me)V^/±D (p, +1 (p ■ g) - Pj-dP ■ «)) f + /rr^) 



;(p) 2 J + 1 



4w(g)* 



(80c) 
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>(P) 



92(p) = h 2 {p)+ l - [ -^V(k) 



2 J (2tt) 3 

M(q) y/J(J+l) 
Q(q) 2J+1 



1 Pj+i(P ■ q) + „ /, : fj-i(p ■ g) ) g2(g) 



2J+1 



2J+ 1 



(Pj+l(p-q)-Pj-l(p-q))9l(q) 



(80d) 



for the four functions 



fci(p) 



w(p) 



(81) 



5i 0) 



(82) 



w(p) 



(83) 



52 b) 



w(p) [fta(p) + 2 ( v /3^ X5 (p) - yZZI X6 (p) 



£1 

4 



(84) 



For J = the Eqs. (|79c|79d|79f1) for X3,X4 and xe are 

absent and all terms in the integrands containing these 
functions vanish in the remaining three coupled integral 
equations. Consequently there exist only the two func- 



tions h2{p) and gi (p) which respect the two coupled in- 
tegral equations (|80c|80d| . 

For mesons of category 3 the BSE becomes a coupled 
system of four integral equations 



xiip) 



Xsip) 



Xiip) 



1 f d 3 q V(k) 
~ 2 J (2tt)3 w 2 (<?) _ i£ 

x \u(q)xi(q) ' 



Pj{p ■ q) 



J 1 



2u(q) 
d 3 q 



J 



2J+1 



X2(q) 



J+l 

2J+1 



XM +M{q) X4 (q) 



2J+12J (2tt)3 W 2 (g ) _ n 



Pj+i{p ■ q) 



u(q) 



j 



j + i 1 



d 3 r 



2J+1 

V(k) 



X2(q) 



2J+12J (2tt) 3 w a(g) _ £ 
1 



Pj-i(p-q) 



J+l 
2J+ 1 

3 



w(g) 



l^ 1 ^ ' 0(g) 

1 /" d 3 g 

27 (2tt)3 ^yrz 

Y Xl(,z) + ^) (/ 



Xa(?) 



2J+ 1 

M(g) 



J+l 
2J+ 1 



J 



2J+ 1 



X2(g) 



J + l 

2J+ 1 



X3(q) +M(q) X i{q) 



XM +M(q)xM 



x 3 (q) +M{q)xM 



(85a) 



(85b) 



(85c) 



(85d) 
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Obviously in the integrands there appear only two linear combinations Xi (l) and the term between the square brackets. 
Thus one can reduce (|85p to a coupled system of two integral equations 



1 f d 3 q M(p)M(q)Pj(p-q)+pq(^ TI Pj +1 (p-q) + ^ I Pj^(p-q) 

u{p)h{p) = - / T^-To^(fc) 



2 J (2^) 3 v ' ufrMq) 

x { k ® + - (86a) 

- 4^)) ff(p) = %) + ^ / (0 ^(fc)Pj(p • <?)<?(<?) (86b) 



for the two functions 



where 



h(p) = 
and 



p 

Hp) 



w(p) 



o;(p) [fe(p) + 2 X i(p)] 



For J = the ansatz for the vertex function contains only 
component 2. Inserting it into the BSE gives equation 
(I85bp which becomes Xz{p) = in this case. Thus a 
meson with quantum numbers 0~~ does not occur in the 
model with instantaneous interaction. 

For mesons of category 4 the BSE becomes 



Xiip) = Xa(p) = 0. 



(89) 



That means that in the model with instantaneous in- 
teraction no mesons of category 4 occur. The quantum 
numbers for J = 2n and J ^ for J = In + 1, re- 
spectively, and also 0~~ of category 3 are called exotic. 



APPENDIX C 

Here we define the rest frame meson wave function and 
project it onto the quark-antiquark component propagat- 
ing forward and backward in time, see [161 ] and references 
cited therein. 

The rest frame wave function is 

iPjm(^P) = 

1 J ^B S(P0 + ^P)X^MP)S(P0 - f , P). (90) 
The quark propagator can be decomposed into 

S(po,p) = S + (po,p) + S-(p ,p) (91) 



S±(po,p) = 



(87) W ith 



Po T w(p) ± ie 



T p = exp 



and the projectors 



-l-py--V P 



(92) 



(93) 



l±7o 



(94) 



Integration over po and using P±7o = ±P± and TJ70 
7 T p yields 



^Pjm(^P) 



t p p+t pX f AMp)T p P-T p 

2u[p) - M 
T p P_T pX p g{p,p)T p P+T p 
2u)(p) + M 



(95) 



Due to Tjj* = T p 1 the Foldy transformed wave function 
tijZM = T^jm(^P)T£ has the form 

Vjm^P) = P+^SmMP- +P-i> PC JM {^P)P+, 

(96) 

with components •4'±jm(^'P) propagating forward and 
backward in time. Notice that for the quark P + and 
P- project out the forward and backward components, 
respectively, but for the antiquark it is in the other way 
around. The general form of ?A^j M (^,p) for a meson of 
given quantum numbers can be derived by introducing 
the respective vertex function into ([9l>]) . This yields 



iP±jm(v>,P) = ■4>±(p)Yjm{p)1\ 



5 — 7: 



Yju{p)l5 
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(97) 



for mesons of category one, 



- \ / - — : 1pl± (p) + \/ o t + \ ^± (P) 



2J+1 



2J+1 



J 



2J+1 



■02±(p) 



{y 7+ i(p)®7} JM 



with 



for mesons of category 2 and 



1 



hi(p) ± 



2 V 2w(p) / 



{yj(p)®7} JM 



(98) 



(99) 



(100) 



for mesons of category 3. 



APPENDIX D 



In this Appendix we show how we numerically solve 
the BSE. We calculate the spectra and wave functions 
with small but finite values of /iiR, specifically for /kir = 
0.01-ycr, 0.005-^ and 0.001y/a. As input we use the re- 
spective results for M(j>) and u>{p) which are obtained 
by an iterative solution of the gap equation. With the 
given small but finite /zir the deviation of a* 2 (/xir) from 
the IR limit /j, 2 = p 2 (0) turns out to be almost linear for 
small values of /zir. Then meson mass fj, in the IR limit 
can be determined by extrapolation to /iir = with high 
confidence. In a similar manner the IR limit of the wave 
functions can be obtained. In Figs. [15] we show as ex- 
ample the determination of the IR limit for the mass of 
the fourth radial excitation of the pseudoscalar meson. 
The forward-propagating component of the wave func- 
tion of the same state is shown in Fig. [TBI together with 
the extrapolated IR limit, which in this plot is hardly 
distinguishable from the result with /iir = 0.001^/a. 



The systems of integral equations ([76]) . ([80]) and 
([86]) are numerically solved by expanding the functions 
h, g (respectively h±, gi, hi, (72) in a finite set ((fi,i = 
1, . . . , N) of basis functions. In the following we show 




FIG. 15: Results for the square of the mass of the fourth radial 
excitation of the pseudoscalar meson in units of a for three 
values 0.01 ^/a, 0.005 v / o : and O.OOIa/ct of the IR regulator pm. 
and extrapolation to the IR limit fim = 0. 



details for mesons of category 1: With 



i = l 

N 



i=l 



(101a) 
(101b) 



one ends up with 



7fJ 



B(n) 



5* 



(102) 
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.... m R =o.oi 




.-. H m = 0.005 




_- H ]R = 0.001 




— t*IR= 






FIG. 16: Component ij)+(jp) of the wave function of the fourth 
radial excitation of the pseudoscalar meson for three small 
values of /iir and its extrapolation to /xir = 0. All quantities 
are given in appropriate units of \fa. 



with symmetric real N x N matrices A hh , A 99 (/i), B hh 
and B^(n) and real N x N matrices B h "(fj,) and B ah (fx) 
with B hg (fi) T = B 9h (fi). The matrix elements are given 
by 



d 3 p 
J2^f 



(pi(p)oj(p)<Pj(p), 



(105) 



d 3 p 
(2tt) 3 



uj(p) 



fj(p), (106) 



where the real symmetric 2N x 2N matrices A(/j,) and 
B(fi) have the block structure 



A{p) = 
B(fi) = 



A 



hh 







B 



hh 



\ j 

( B hh B hg (ii) 
{B" h (p) B"»(jjl) 



(103) 

(104) 
I 



1 f d 3 pd 3 q 



^{p)V{k)Pj{p-q) V] {q)^ (107) 



(108) 



' M(p)M{q)Pj(p -q)+pq (^Pj +1 ( P ■ q) + j^jPj-iiP ■ g) 



w(p)w(g) 



4w(p)w(g) 



^j(p ■ g) 



^•(g). (109) 



One can drop or add terms of order 0(/LtiR,) which leads 
to different results at finite /zir but leaves the IR limit 
invariant (see the discussion on the IR limit of the BSE 
above). For instance one can drop the /i-dependent terms 
in A aa {n) and B S9 (/x) or replace ([108]) by 

B??(P) = B#G*) = § / J^i(pVi(p)- (HO) 

We have checked that these modifications indeed give — 
up to insignificant numerical errors — the same extrap- 
olated IR limits. All results presented in this paper have 
been calculated with the exact expressions (|105l) - (|109p 



for finite /iiR. By introducing a multiplicative param- 
eter A(/i) on its right hand side, Eq. (|102[) becomes a 
generalized eigenvalue problem which can be solved with 
standard linear algebra methods yielding 2N eigenvalues 
Aj(/x). The masses /Ltj (i — 1,...,K) of the K lowest 
lying mesons with given spin are then found by solv- 
ing Xi(fii) = 1 for the lowest K eigenvalues. The corre- 
sponding eigenvectors give the coefficients in the expan- 
sion (jlOip in each case. Finally the functions h and g 
have to be normalized according to (|4"2"|) . The spectra 
and wave functions of mesons of categories 2 and 3 are 
found in an analogous manner. As basis functions we 
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take 



¥>i(p) = P L exp(-«ip 2 ) 



(111) 



where L = J for mesons of category 1 and 3 and 
L = | J — 1| for mesons of category 2. With properly 
chosen parametrs a number of N — 20 basis func- 
tions is sufficient for a satisfactory numerical accuracy of 
masses and wave functions of all states considered in this 
work. 
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